% Plot confocal scan image

function [h,ax] = plot_confocal(image_file,fontsize)

[IMG, ~] = ReadImageFile(image_file);
ax = axes;
xvals = linspace(IMG.minVx,IMG.maxVx,IMG.NX);
yvals = linspace(IMG.minVy,IMG.maxVy,IMG.NY);
[xvals,yvals] = meshgrid(xvals,yvals);
h =pcolor(xvals,yvals,IMG.Data);
daspect([1 1 1]);
colormap(ax,pink);
shading interp
b = colorbar;
ylabel(b,'Count rate (Mcps)','fontsize',fontsize)
set(ax,'fontsize',fontsize)
xlabel('X(um)','fontsize',fontsize)
ylabel('Y(um)','fontsize',fontsize)